Accepted for publication by the Astrophysical Journal 

Preprint typeset using L^T^X style cmulatcapj v. 5/2/11 



THE INITIAL MASS FUNCTION OF THE ORION NEBULA CLUSTER ACROSS THE H-BURNING LIMIT 

N. Da Rio 

European Space Agency, Keplerlaan 1, 2200 AG Noordwijk, The Netherlands 



M. ROBBERTO 

Space Telescope Science Institute, 3700 San Martin Dr., Baltimore MD, 21218, USA 



L. A. Hillenbrand 

California Institute of Technology, 1200 East California Boulervard, 91125 Pasadena, CA, USA 



T. Henning 

Max-Planck-Institut fur Astronomic, Konigstuhl 17, 69117 Heidelberg, Germany 



K. G. STASSUN 2 ' 3 

Vanderbilt Univ., Dept. of Physics & Astronomy 6301 Stevenson Center Ln., Nashville, TN 37235, USA 
Accepted for publication by the Astrophysical Journal 

ABSTRACT 

We present a new census of the Orion Nebula Cluster (ONC) over a large field of view (> 30' x 
30'), significantly increasing the known population of stellar and substellar cluster members with 
precisely determined properties. We develop and exploit a technique to determine stellar effective 
temperatures from optical colors, nearly doubling the previously available number of objects with 
effective temperature determinations in this benchmark cluster. Our technique utilizes colors from 
deep photometry in the /-band and in two medium-band filters at A ~ 753 and 770 nm, which 
accurately measure the depth of a molecular feature present in the spectra of cool stars. From these 
colors we can derive effective temperatures with a precision corresponding to better than one-half 
spectral subtype, and importantly this precision is independent of the extinction to the individual 
stars. Also, because this technique utilizes only photometry redward of 750 nm, the results are only 
mildly sensitive to optical veiling produced by accretion. Completing our census with previously 
available data, we place some 1750 sources in the Hertzsprung-Russel diagram and assign masses and 
ages down to 0.02 solar masses. At faint luminosities, we detect a large population of background 
sources which is easily separated in our photometry from the bona fide cluster members. The resulting 
initial mass function of the cluster has good completeness well into the substellar mass range, and we 
find that it declines steeply with decreasing mass. This suggests a deficiency of newly formed brown 
dwarfs in the cluster compared to the Galactic disk population. 

Subject headings: stars: formation, pre-main sequence, mass function, brown dwarfs, stellar clusters: 
individual (Orion Nebula Cluster) 



1. INTRODUCTION 

Understanding the initial mass function (IMF) is one of 
the most important problems of stellar astrophysics. The 
IMF, together with the star formation history, dictates 
origin, evolution, and fate of the stellar populations, from 
individual clusters up to entire galaxies. The distribu- 
tion of stellar masses has been studied in depth for more 
than half a cen tury, starting with the pioneering work of 
Salpeter (1955); the major open question is whether the 
IMF is universal or if it dep ends on the in i tial conditions 
of sta r formation (see, e.g., Kroupa| |2002 ; Bastian et al. 
20101. Whereas it is commonly accepted that the IMF 



seems well reproduced by a power law for masses greater 
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than several tenths of a solar mass, the shape and uni- 
versality of the IMF in t he substellar mass regime is still 
under investigation (e.g., Wang et al. 2011). 

Very young clusters (few Myr old) m star-forming re- 
gions provide a unique tool to investigate the IMF across 
the entire mass spectrum, for a number of reasons. First, 
they are usually young enough that neither dynamical 
processes nor stellar evolution have altered the mass dis- 
tribution; therefore the measured distribution of stellar 
masses coincides with the IMF. Moreover, young low- 
mass stars and brown dwarfs (BDs) are in their bright- 
est evolutionary stage when they contract towards the 
main sequence, therefore these objects are more easily 
detected and characterized when they are young. 

Among the nearby Galactic star-forming regions, the 
Orion Nebula Cluster (ONC) is an ideal site for the study 
of star formation in particular for low-mass stars and 
BDs. This cluster counts a few thousand members, 1- 
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3 Myr old ( |Hillenbrand|[l997l |DaTRio et al |[20lo| , span- 
ning the entire mass spectrum [M < 5U Mq). Several 
studies have been conducted i n the past decade to mea- 
sure the IMF in the ONC (e.g.,|ffillenbrand|1997||Hillen-| 
brand fc Carpenter|2000l|Muench et al.|2000J|2002||Lucas1 
fc Roche|2000| [Lucas et al.|2005HSlesnick et al.|2004||Da| 
Rio et al.||2010 ); they generally find (similar to the field 
star population) a Salpeter-like slope above 1 M , which 
flattens to a broad peak at 0.2-0.3 M (though the shape 
and position of the IMF peak is highly model dependent, 
Da Rio et al~|2010[ ), and the mass distribution likely de- 
creases in the substellar mass range. 

Determining masses in a young region such as the ONC 
is not without difficulties: besides the strong nebular 
emission, another major impediment is caused by dif- 
ferential reddening, which has strongly limited the abil- 
ity to derive the stellar parameters of individual sources 
based on photometry alone. To overcome this shortcom 
ing, spec trosco pic surveys h ave been ca rried out flHil 
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Riddick et al.||2007[ | Weights et al.||2009l ), but they 
are either limited to a fraction of the members, or to 
the very central part of the region, the Trapezium clus- 
ter. Using near-infrared (NIR) photometry it is possi- 
ble to asses stellar masses down to the planetary masses 
(M < 13 Mj), for example by de-reddening the measured 
color- magnitude diagrams (CMDs) on one isochrone or 
even to deriving directl y the IMF from the NIR lumi- 
nosity functions (L Fs) (Hill enbrand fc Carpenter||2000[ 
Muench et al.|[2002 ). This second approach, however, is 
not very accurate: the ONC shows a significant luminos- 
ity spread, sometimes interpreted as a real age spread 
(Reggiani et al. 2011), and stars of different age follow 
different mass-lummosity relation. Moreover the NIR ex 



cess origina ting in the inner circumstellar disks ( |Meyer 
et al . |[l997 > alters the observed fluxes 



Clearly, in order to improve our knowledge of the ONC 
IMF, a precise and systematic characterization of the 
stellar parameters of individual members is neede d. To 



this purpose, the Orion HST Treasury Program (Rob- 
berto et al.|2005 ) has produced a high spatial resolution 



photometric survey of the ONC with three instruments 
onboard the Hubble Space Telescope, from U— band to 
the NIR, over a large field of view (FOV) , with a sensi- 
tivity well into the BD mass range. To fully exploit this 
exceptional dataset, an accurate estimate of the effective 
temperature (T c g) is needed to derive Ay from the ob- 
served colors and t herefore the stellar lumino sities. In 
our previous works (|Da Rio et al.|[2009b 2010 hereafter 
Paper I and Paper 11) we have collected spectral types 
from the literature and complemented them with new op- 
tical spectroscopy. Moreover, we have presented a new 
observational technique, based on optical medium-band 
photometry at 6200A — a wavelength where the spec- 
tra of cool stars show a deep, T e ff-dependent TiO ab- 
sorption feature — to derive the spectral types of M-type 
sources. This allowed us to obtain the stellar parameters 
of ~ 1000 members and derive the most complete HRD 
of the ONC, down to the hydrogen burning limit. 

In the present paper we extend our investigation to 
lower masses, well into the BD regime, and with higher 
completeness. To this purpose, we take advantage of 
the same observational strategy presented in Paper I, 



using medium-band photometry to derive spectral types 
of cool sources, from which T e g and Ay are derived, and 
place individual sources in the HRD. Here, instead of the 
6200A filter used in Paper I, we select two bands at longer 
optical wavelengths, in order to increase our sensitivity 
for cooler sources such as BDs and further reduce the 
influence of extinction. 

In Section[2]we present our new observations, the data 
reduction and calibration. In Section [3] we define two 
spectro-photometric indices based on our medium-band 
photometry. Using available stellar parameters for a 
small sample of ONC very-low mass stars and BDs, we 
define empirical transformations to convert these indices 
into T c ff and Ay. We also show that this technique is 
not significantly affected by optical excesses from mass 
accretion, typical of young stars and BDs. In Section [4] 
we derive the new HRD for the ONC, which now includes 
~ 1800 sources, reaching masses as low as 0.02 M . We 
study the completeness as a function of T e g and bolomet- 
ric luminosity (ibol) accounting for photometric detec- 
tion as well as differential reddening; we detect a popu- 
lation of candidate background contaminants, which ap- 
pears well separated from the ONC members in the HRD, 
at lower luminosity. After excluding the contaminants, 
in Section [5] we derive and discuss the IMF of the ONC. 

2. THE DATA 
2.1. Observations 

Observations were carried out with the Wide Field Im- 
ager (WFI), a focal reducer- type camera mounted at the 
Cassegrain focus of the 2.2-m MPG/ESO telescope at La 
Silla. The field of view (FOV) of the camera is - 32' x 33' , 
allowing us to cover the entire ONC with one pointing. 

This optical imager offers a particularly large set of 
broad-, medium-, and narrow-band filters along the en- 
tire wavelength range of its CCD sensitivity (~3000 to 
lOOOOl). We selected two WFI medium band filters 
which provide the best sampling of a deep absorption 
band at A - 77001. These are the MBj}753/18-ES0848 
and the MB|j770/19-ESO849 filters, which hereafter will 
be simply referred to as 755 and 770. Figure [T] shows 
the transmission pro files of these fil t ers ov erlayedto syn- 
thetic spectra from Allard et al. (2000) of decreasing 
temperature. The 753 band is centered on the photo- 
spheric continuum, while the 770 is completely covering 
the TiO absorption band whose strength increases with 
decreasing T e g. We also used a broad /-band filter, the 
BBjjI/203-ESO879, to sample the stellar flux at redder 
wavelengths, in order to evaluate the photospheric red- 
dening from dust extinction. This I-band filter is also 
used in our previous photometric studies of the ONC 
(Paper I and II). Table [l] summarizes the characteristics 
of the 3 filters. 

The WFI imaging was carried out during 5 nights of 
December 2010. On each night, observations were per- 
formed in all the three bands together so as to mini- 
mize the impact of any longer timescale stellar variabil- 
ity. A standard dithering pattern has been applied, in 
order to cover the gaps between the 8 CCD chips of WFI 
and allow for cosmic ray and bad pixel removal. We 
did not carry out short exposures to avoid saturation for 
the bright sources, since the bright end of the ONC has 
been studied well in detail by us in our previous works. 
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Fig. 1. — The band profiles of the 3 WFI filters used in this 
work (from left to right: 753, 770 and I — band) overlaid to spectra 
of very low- mass stars (4000 K t o 3000 K) and BDs (3000 K to 
2500 K) from |Allard et al-lpOOO) . 



TABLE 1 

Characteristics of the filters used for our WFI 

IMAGING 



band 


A central 

(A) 


FWHM 

(A) 


E.W. a 

(A) 


Vega zero-point 
(erg s -1 cm -2 A -1 ) 


753 


7536.2 


184.1 


191.4 


1.368-10" 9 


770 


7707.9 


195.2 


201.3 


1.275-10- 9 


I 


8620.5 


1354.9 


1431.2 


9.500-10- 10 



Equivalent Width 



TABLE 2 

WFI OBSERVATIONS 



band 


exposure 

w 


time 


n. of exposures 


total exposure time 
(s) 


753 




360 


23 


9,360 


770 




500 


28 


14,000 


I 




280 


12 


3,360 



Dark, bias and flat field exposures were obtained before 
and after each night. We did not observe standard fields 
for photometric calibration, given the particular filters 
used in this program and the practicality of performing 
relative photometric calibration exploiting our large pre- 
existing data set on Orion. 

2.2. Data reduction 
Images were processed using the E SO/MVM (Multi 



Vision Model) vers. 1.3.5 package ( Vandame ||2004[ ). 
This software automatically performs all common image 
reduction steps (e.g., correction for bias and flat fields), 
and the absolute astrometric calibration of the individ- 
ual exposures to be merged. This is accomplished by 
cross-correlating the images with a reference astromet- 



ric catalog; we utilized to this purpose our previous WFI 
photometric catalog presented in Paper I. The final prod- 
uct of our image reduction is 3 images, one for each filter, 
obtained by merging all the individual exposures. Table 
[2] summarizes the total exposures times; these about 1 
hour for the broad band, and about 2.5 and 4 hours for 
the two medium bands, 

Photometry, both aperture and point-spread-function 
(PSF), was p erformed using the Daophot II package 
(|Stetson||1987|). We computed the PSF individually for 



every image, starting from a coadded sample of few hun- 
dreds bright and unsaturated sources, and then refining 
these with an iterative sigma-clipping method to reject 
PSF stars with poor \ 2 - Photometry was computed on 
all the sources brighter than 3<r above the local sky back- 
ground noise. In order to eliminate the contamination 
from spurious detections — a potentially significant frac- 
tion of the total detections, given the strong variability 
of the nebular background, and the low threshold for 
source detection — we rejected sources not present in any 
of our existing deep photometric catalogs of the ONC, 
both grou nd based (WFI UBV I, Paper I; CTIO/ISPI 
JHK, see |Robberto et~aT1|2010| ) or fro m the HST pho- 
tome try (ACS BV1Z, WFPC2 UBI, see|Robberto"et~aT 



2005]) . In particular, the significant depth of the infrared 
and HST data, over a FOV equal or larger than that of 
the observations described here, guarantees that we are 
not excluding any real sources in this study. 

The photometry has been calibrated to the VegaMag 
system using our previous WFI photometric catalog. For 
the /-band, we simply determined the photometric off- 
set with respect to our previous photometry in the same 
filter. For the two medium bands, we used synthet ic pho- 
tometry. We considered the atmosphere models of |Allard] 



et al. ( 2000 ) as empirically calibrated by us in Paper II, 
and computed the two medium-band colors (MB — I) as 
a function of T g. Then we considered the actual sources 
detected in all 3 bands and for which we have T c g and 
Ay from Paper II, and plotted the measured, extinction- 
corrected (MBi nstrumon t a i-/ ca /i 6 ,. Qted ) colors against T cS . 
The systematic difference between these colors (consid- 
ered individually) and the synthetic ones is the offset 
needed to calibrate the medium bands. Specifically, we 
limited to the stars with T c g > 4000 K, range where the 
possible inaccuracies of the atmosphere models are very 
modest. 

2.3. Completeness 

We derive the completeness functions for our WFI pho- 
tometry with an artificial star test, performed separately 
for each WFI filter. Artificial stars of different magni- 
tude, are added on the reduced WFI images; then, with 
the same technique used to detect and extract photome- 
try on the real stars, these artificial stars are recovered or 
missed depending on their flux, local sky noise, crowding, 
etc. The fraction of recovered stars against the total, as 
a function of magnitude, provides the completeness. 

In our ONC imaging, the completeness is spatially vari- 
able, being strongly affected by the nebular background, 
which is brighter at the center where the stellar density 
is higher. This correlation, if not accounted for, may 
bias the derived completeness. Therefore, we first derive 
the spatial, projected, stellar density distribution of the 
ONC. This is obtained by merging all catalogs at our 
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Fig. 2. — Photometric completeness functions for our three filters. 
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Fig. 3. — Photometric errors as a function of magnitude, for the 
3 photometric bands 

disposal, in particular the deep ISPI JHK photometry 
and the HST/ACS data. Then, for every magnitude bin, 
the artificial stars are generated with random positions 
drawn from the 2 dimensional density distribution. 

The derived completeness functions are shown in Fig- 
ure [2j In all cases, these curves show a characteristic 
trend: besides the typical steep cutoff at faint luminosi- 
ties which traces the detection limit, the completeness 
shows a slow decrease below 1 at bright luminosities. 
This is due to the luminous nebular emission and sig- 
nificant crowding of the central part of the ONC, where 
the detection limit is significantly poorer. In the outer 
regions of the ONC, we are able to detect sources all the 
way down to»n~ 22.5 mag in / band (see Figure[4]), and 
tom~ 23.5 mag in 753 and 770. 

2.4. Photometry results 

The final catalog includes 2474 sources with photome- 
try in all 3 bands. The photometric errors as a function 
of magnitude are shown in Figure [3] Figure [4] shows our 
new CMD, highlighting the significant improvement in 
depth with respect to our previous work (Paper II). In 
particular, our photometry covers a very large range of 
stellar luminosities, spanning ~ 9 mag in /-band, and 
the new observations extend ~5 mag deeper than in the 




(753-1) 

Fig. 4. — Open circles: the derived I vs 753 — / CMD including 
only the sources measured in all the 3 bands (i.e. also in 770). The 
filled blue circles indicate the sources included in our previous study 
of the stellar population of the ONC, based on WFI photometry 
and spectroscopy down to the H-burning limit (Paper II). The 
average photometric errors as a function of J-band are indicated 
to the right. The red line represents a 1 Myr |Baraffe et al.|^1998) 
isochrone; the square, circle and triangle symbols along the line 
mark masses of 0.4 M Q , 0.08 M and 0.02 M n. The arrow denotes 
the reddening vector for Ay = 2 assuming the |CardeIli et al.| | fl989l l 
reddening law for Ry =3.1. 

previous work, reaching I ~ 22, which corresponds to 
young sources with masses < 0.02 M at 1 Myr (see Fig- 
ure 13 . In Figure E] we also observe among the faintest 
sources a concentration of objects with I > 18 which, as 
we show below, comprise a distant, background stellar 
population that is detectable despite the high extinction 
through the molecular cloud because of the depth of our 
observations. 

In Figure[5]we present the (753—/) vs (753—770) color- 
color diagram, which forms the basis of our analysis in all 
that follows. The most salient feature of this diagram is 
the inflection point near (0.05,0.25), from which the stars 
extend in two different directions. By comparison with 
the orientations of the synthetic isochrone and the ex- 
tinction vector, we see that this inflection in the diagram 
essentially separates the stars into a group (extending to 
the left) that stretches mostly parallel to the isochrone 
and perpendicular to the reddening vector, and a sec- 
ond group (extending to the upper right) that stretches 
parallel to both the isochrone and the reddening vector. 
The inflection in the diagram occurs at T s w 4000 K. 

For stars with T eS < 4000 K (i.e. spectral types >M0), 
the (753 — 770) color is both highly sensitive to T e g, be- 
coming rapidly bluer (more negative) as T e ff decreases 
(see also Figure [I]), and largely insensitive to Ay be- 
cause of the small wavelength separation between the 
753 and 770 bands. In our analysis below we exploit 
this property of the diagram to break the degeneracy be- 
tween T G g and Ay for these very cool stars, which has 
previously been one of the most difficult challenges in the 
derivation of stellar parameters fr om optical broadb and 
photometry for M-stars and BDs ( Hillenbrand||1997[ Pa- 
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Fig. 5. — (753 — 7) vs. (753 — 770) color-color diagram from our photometry. The color and size of the dots indicates the I-band luminosity 
of each source, as ind icated. The black arr ow indicates the direction of the reddening vector for 2 mag of extinction in V— band, assuming 
the reddening law of |Cardelli et ah] | |1989[ l for Ry = 3.1. For refere nce, the thick red lin e is the empirical isochrone derived in Paper II 
converted in this color-color plane using the atmosphere models of |Allard et al.| (|2000). For guidance, the origin (0,0) corresponds to 
10,000 K, and the turn-over in (753 — 770) corresponds to T e gf ~ 400U K, i.e. between K and M spectral classes. 



per II) . We note here that the synthetic isochrone clearly 
does not correctly reproduce the slight upward curva- 
ture of the stars as one moves to the coolest T c g, but 
because these cool stars nonetheless extend perpendicu- 
lar to the reddening vector, we demonstrate below that 
the systematic error in the (753 — I) colors predicted by 
the synthetic isochrones can be readily calibrated out for 
T e s < 4000 K. We note also that the effects of contam- 
ination by faint background sources is mitigated in this 
region of the diagram, particularly for the coolest ob- 
jects of interest (T c g < 2900 K, i.e., in the BDs regime), 
since background BDs are too intrinsically faint to pen- 
etrate the very high Ay of the cloud behind the ONC. 
As we show in Section [4j the relatively few remaining 
background contaminants in this region of the diagram 
(mostly red giants) can be filtered out on the basis of 
their luminosity. 

These desirable features of the (753 — I) vs (753 — 
770) color-color diagram break down for warmer stars 
with T e g > 5000 K. For these warmer stars the colors 
are largely degenerate in T e g and moreover the isochrone 
and the reddening vector are parallel. For these stars 
we must therefore rely on previous observations by us 
and from the literature to determine T e g- and Ay. In 
addition, there is significant background contamination 
in this region of the diagram; as we show below, while the 
brightest of these stars are bona fide ONC members, the 
faintest are mostly background stars. Finally, we note 
that the apparent lack of low- Ay stars for T e g > 4000 K 
is the result of these stars being saturated in our deep 



images. Such stars are indeed present in the ONC, and 
we include them in our analysis from our own previous 
work (Paper I and II). 

3. ANALYSIS 

In order to derive the stellar parameters of very low- 
mass stars (VLMSs) and BDs in the ONC, we focus our 
analysis on the left-hand part of the color-color diagram 
of Figure |5j where the (753 — 770) color anti-correlates 
with T cff . 

Our basic methodology is to project each source from 
its observed position in the color-color diagram back 
along the reddening vector to the theoretical isochrone. 
The choice of theoretical isochrone is largely arbitrary 
because none of the available isochrones exactly repro- 
duce the observed distribution of the stars in the dia- 
gram, and consequently we must empirically recalibrate 
the isochrone based on the observed distribution. Our as- 
sumption is that, whichever fiducial isochrone is adopted, 
the non-reddened stars should lie on the isochrone, and 
the reddened stars displaced away from it along the red- 
dening vector. For example, the isochrone shown in Fig- 
ure [5]computed using the synthetic spectra of | Allard~et] 
|aL| |2000[ ) predicts a much bluer (753 — 7) index than 
our data for low (753 — 770). If the model were correct, 
this would imply an unrealistic dependence of the aver- 
age Ay on T e g , in the sense that the minimum measured 
extinction (i.e. the distance, along the reddening vector, 
between the isochrone and the stars) increases monoton- 
ically with decreasing T e g. 

Different synthetic isochrones can lead to more realistic 
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Fig. 6. — The color-color diagram of Figure [5| with different 
synthetic isochrone s overplotted. In part icular: a piess et al.||2000]l 
2 Myr isochr one; alBaraffe et al.| 11998} 1 Myr isochrone; alsiessl 
let all pOOOl zero-age main sequence. Thick lines are the result 
from synthetic pho tometry on the AMES-MT grid of spectra from 
|Allard et al.| ( |2000|) , while d ashed line are computed using the BT- 
Settle grid from |Allard et al.||2bl0} . The thick red line is computed 
using the empirical calibration ot the AMES-MT models defined 
in our Paper II to match the broad band BVI colors of the ONC. 
The empirical line is actually closer to the predicted ZAMS relation 
than the predicted 1-2 Myr isochrones, illustrating the current lack 
of agreement between evolutionary and atmosphere models and 
young star photometry (see text). 
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Fig. 7. — Conceptual scheme of the new coordinate system we 
introduce to measure T e g and Ay of the stars from our color- 
color diagram: the T e ff dependent [770 index], and the reddening 
dependent Ay index. See text. 

predictions of colors. Figure [6] shows our color-color di- 
agram with a number of synthetic is ochrones, computed 
using both the AMES-MT spectra of |Allard et"aI] ( |20Q0[ ) 
and the most recent grid — the BT-Settle models— from 
Allard et al. (2010). We also show an i sochrone (red 



line] derived using the Allard et al. ( 2000[ ) atmospheres 
as recalibrated by us in Paper 11. That calibration, ap- 
propriate for the ONC, involved an empirical constraint 
on the stellar log g to produce a synthetic isochrone that 
matched the observed broadband colors of the ONC in 
BVI. It is clear that all synthetic isochrones, including 
that which we used in Paper II, are incompatible with the 
data in this new medium-band color-color plane. AMES- 
MT consistently predicts a systematically lower (753 — /) 
at low T e g, and BT-Settle predicts either a turn-over in 
(753 — 770) for low log g (young isochrones) or (753 — I) 



colors that are too red for high log g (the ZAMS model) . 
Evidently, current synthetic models, even those that have 
been empirically calibrated by us to match the broad- 
band colors of young low-mass stars, are unable to cor- 
rectly reproduce the molecular feature at A ~ 7700A 
and/or the photospheric continuum level at A ~ 7500A 
which form the basis of our new medium-band color-color 
diagram technique. As a consequence, we must empir- 
ically define an isochrone representing both the photo- 
spheric colors in our bands and the T e g scale along it. 

To proceed, we introduce a new set of coordinates 
in the color-color plane of Figure [3J We define the 
[Ay index] as the distance that a given star is displaced, 
along the reddening vector, from the x— axis of the color- 
color diagram, measured in units of Ay : 



[Ay index] 



m 7 53 



A 7 S3 



Al 

Av 



(1) 



where A 753 /A v = 0.678, A 770 /A v = 0.652 and A T /A V = 
0.533 (values computed from the |Cardelli et al" (1989) 
reddening law for Ry = 3.1). We choose the redden- 
ing parameter Ry = 3.1 instea d of a higher value (e.g, 
5.5, Costero & Peimbert 1970) because, as we found in 
Paper 11, this value better explains the observed broad- 
band colors of the ONC members. We similarly define 
the [770 index] as the projection of a given star, along 
the direction of the reddening vector, onto the x-axis [i.e. 



(753-7) = 0]: 
[770 index] =771753 
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By definition, the [770 index] is reddening indepen- 
dent, and so it depends only on T e g. Figure [7] shows a 
schematic representation of the two indices in the color- 
color diagram. 

3.1. Derivation of T c g 

Because the currently available synthetic spectra are 
unable to accurately predict the stellar colors in our 
medium-band color-color diagram, we must derive an 
empirical transformation between the [770 index] and 
T ff. To this end we utilize the sample of ONC stars 
with directly determined (i.e. spectroscopic) spectral 
types. With the appropriate transformation between the 
[770 index] and T e ff, we will be able to reliably deter- 
mine T e ff of all the other members from their measured 
[770 index] . We have assembled all available optical spec- 
tral types from our Paper II (for stellar-mass objects) 
and NIR spectral typ es ( for substellar-mass ob jects) from 
Riddick et aL] ( |2007[ i and |Slesnick et al.| ( |2004[ ). The spec- 
tral typ es were converted to 'Jeff using the temperature 
scale of Luhman et al. (2000), as we did in Paper II. 

Figure IS] shows the resulting relationship between our 
measureci[770 index] and the spectroscopically deter- 
mined T c ff. As expected, for late spectral subtypes the 
[770 index] becomes bluer (i.e. becomes more negative) 
with T c ff. This approximately linear dependence holds 
from T cff ~ 3500 K (~M2.5) down to the coldest avail- 
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Fig. 8. — Derived relation between the measured, reddening in- 
dependent, photometric [770 index], and spectroscopically deter- 
mined stellar temperature. The values of T e g for the stellar mass 
regime (T e fj- > 3000 K) are from the spectral types collected in our 
Paper II. For these sources, we identify a sample of candidate weak 
accretors, showing low Ha emission (E.W.< 10A) and no signifi- 
cant veiling (log L accI /Ltot < —1-5), in order to more cleanly de- 
termine the true relationship between the photometric [77 index] 
and T.ff . Sp ectr oscopic T„fr for substellar masses are from |Riddick| 
|et al.|{2007] and |Slesnick et al.| |2004] |. See text for details. 

able spectral types for BDs, and down to our photometric 
lower limit ([770 index] ~ 0.8). 

Perhaps not surprisingly, we found that the scatter 
in the [770 index]— T e g relation could be significantly 
reduced if we consider only sources that do not show 
strong spectroscopic evidence of accretion. We selected 
the "non-accretors" on the basis of low Ha emission 
(E.W.< lOA) from Paper I and low optical veiling (low 
B— band excess emission) from Paper II. Figure [8] shows 
that these non-accretors form a relatively tight sequence, 
while the rest of the stars (the accretors) are much more 
scattered in the diagram. The accretors are also system- 
atically shifted to higher T e g and/or lower [770 index], 
indicating that accretion alters the position of the sources 
in the diagram, and indeed below we exploit this feature 
to quantify the degree to which accretion may affect our 
determination of T c g. 

Using only the sample of non-accretors, we derive an 
empirical linear relation between T e g and the [770 index] : 



T cff = 3620 K + [770 index] • 1150 K 



(3) 



which is valid for [770 index] < —0.1. The standard de- 
viation of the data relative to this relation is ~ 100 K, 
roughly corresponding to one-half of a spectral subtype. 
This is comparable to the overall uncertainty in our spec- 
troscopically determined T e s from Paper II, suggesting 
that this scatter is principally due to uncertainty in the 
(spectroscopic) spectral types, and therefore any addi- 
tional uncertainty introduced through the [770 index]- 
T e ff transformation is negligible. 

3.2. Influence of accretion excess on the T e g estimation. 

The analysis methodology utilizing our medium-band 
color-color diagram is predicated on the assumption that 
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Fig. 9. — Same as Figure [8] showing the displacement in the 
reddening-independent [770 index] obtained by adding an amount 
of accretion luminosity to the star. The thick line is the isochronc 
from Paper II, the same model shown also in Figure[5] The dashed 
line are the same model, after the addition of a given amount of 
accretion luminosity to the each synthetic spectrum, for different 
values of log L accI /L to t ■ Note that the fact that the red dots are 
low accretion sources, whereas the model requires a large value of 
iaccr to predict their colors, implies that the base isochrone needs 
adjustments (see text). 

the stellar colors in our photometric bands depend only 
on T e ff and Ay, permitting us to establish an empirical 
relationship between the our [770 index] and T c q (Equa- 
tion I3J. However, the colors of young stars can be af- 
fectea by accretion, whereby infalling matter from the 
circumstellar disk onto the stellar surface leads to a flux 
excess, strongest in the UV and in specific emission lines, 
but also to a smaller extent in the continuum at optical 
wavelengths. The latter is typically characterized by a 
heated photosphere covering ~1% of the stellar surface 
and with T R ff ~ 6000 — 8000 K, th erefore peaking in the 



B-band (Calvet & Gullbring 1998). Since our photome- 
try uses filters at the red end of the optical spectrum, and 
there are no emission lines within the wavelength range 
covered by our filters, we expect that accretion should 
induce only a modest alteration of our measured colors. 

We investigate this hypothesis in a quantitative way, 
calculating the shift in [770 index], at constant T e ff, ob- 
tained by contaminating the photospheric fluxes with 
varying degrees of accretion luminosity. As in Paper II, 
where we performed a similar modeling to derive the 
change in the BVI colors due to mass accretion, we first 
model a typical accretion spectrum, then we add it with 
increasing relative flux fraction on synthetic stellar spec- 
tra of different T c ff . Finally, by means of synthetic pho- 
tometry as above, we compute the resulting changes in 
the [770 index] and in turn the error induced in the in- 
ferred T c g. Specifically, we simulated a typical accretion 
spectrum using the Cloudy photoionization code. This 
is similar to that we used in Paper II, but has been fur- 
ther refined to provide a more realistic SED in the ultra- 
violet range (Manara et al. 2011, in prep.). We use the 
synthetic spectra that we used in Paper II-namely the 
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Fig. 10. — Computation of the offset on the derived T c ff, as a 
function of temperature, produced by neglecting accretion veiling. 
The solid lines represent the result for different accretion luminosi- 
ties. The dotted lines represent the uncertainty in T e ff associated 
to un uncertainty in the spectral type of 1 subtype, half subtype 
and 1/4 of a subtype. 



em pirical ly calibrated AMES-MT model from Allard et 



al. ( 2000 ) — and for every T c g we add a component of our 



accretion spectrum to the photospheric spectra, varying 
the ratio L accr /L star . 

The result is shown in Figure [9] from which two main 
results are evident: First, contaminating the stellar flux 
with accretion excess shifts the [770 index] to higher val- 
ues, and the effect is more prominent for cooler T e g. 
This is intuitive considering that this index traces the 
ratio between the flux outside (A ~ 7530A) and inside 
(A ~ 7700A) the TiO molecular feature. The accretion- 
induced continuum excess partially "fills in" the TiO ab- 
sorption feature, leading to higher (less negative) values 
of [770 index). Second, for typical values of accretion 
luminosities in a PMS cluster (£ acC r ^ -^star) the shift in 
the [770 index] is negligible. 

From Figure [9] it also appears that the synthetic 
isochrone does not quite follow our data. This is not sur- 
prising, given that, as discussed in Section [2] this model 
provides a somewhat inaccurate estimate of the stellar 
fluxes in our medium bands (see Figure [6]). Therefore, 
in order to obtain a correct, quantitative assessment of 
the accretion-induced shift in the [770 index], we must 
recalibrate the results shown in Figure |9| We proceed as 
follows. The systematic error in the depth of the TiO 
feature as a function of T e g can be treated as a simple 
systematic error in the T e g scale. For example, for the 
lowest T off modeled in Figure M (T cS ~ 2950 K), the 
models predict a depth of the TiO feature correspond- 
ing to [770 index]= —0.8 mag, and by adding 10% of 
accretion luminosity one obtains a shift of the index of 
~ 0.06 mag. In fact, from Figure [8] and Equation [3j 
[770 index]= —0.8 mag corresponds to a lower T e g, about 
2700 K. Therefore we can associate this T e g with a shift 
of ~ 0.06 mag. Then, we readjust the ratio £ accr /L sta r, 



using the updated T e g, according to L star oc T* s . 

Finally, we can recast the result from a more observa- 
tional perspective, in which we measure [770 index] for a 
given star and we do not know the true value of L acC r for 
that star. Deriving T c g via Equation [3j which is defined 
for L accr = 0, we will overestimate the true T e g because 
accretion will cause the TiO feature to be filled-in and 
con sequ ently we will measure a higher [770 — index]. Fig- 
ure [lO] shows the resulting systematic overestimation of 
the derived T e g for various L accl . We find again that the 
lower the stellar T e g, the larger the error in the T e g es- 
timate due to accretion. However, if one excludes very 
strong accretors (£ accr ~ Lstar) > the error in T e g is very 
small. For instance, for i accr = 10%L s t ar , the T e g error 
is always smaller than half spectral subtype, and always 
< 40 K above the H-burning limit (T cff - 3000 K). The 
T e g difference associated to an error of 1 spectral type 
was computed as the difference between the temperature 
of two contiguous spectral subtypes. Similarly, the T e g 
differences corresponding to 1/2 and 1/4 spectral sub- 
types (dotted lines in Figure [l0|) is equal to that for 1 
spectral type, scaled down by a factor of 2 and 4. 

In conclusion, the presence of optical excess due to 
accretion has a negligible effect on the stellar parameters 
derived from our medium-band survey. Thus, we will 
neglect accretion in the following analysis, and proceed 
to derive T Q g based on the [770 index] from Equation [3j 

3.3. Accuracy of T c g from previous works 

Despite having shown that accretion excess not not 
alter significantly the measured [770 index], Figure [8] 
showed that accreting ONC sources are displaced in the 
T c g vs. [770 index] plane. The displacement is such that 
the accretors lie above the non-accretors, meaning that 
the spectroscopic T c g plotted on the y-axis of Figure 
have been on average overestimated for accreting ON 
sources. 

This overestimation of spectroscopically determined 
T e g for accreting sources is likely the result of incom- 
plete removal of the optical veiling in the spectra of these 
sources. This excess, that we demonstrated above to be 
minimized at 7700A(alth ough it may still be present in 



extreme cases, see, e.g., Fischer et al. 2011), is likely 



larger at bluer wavelengths, where other molecular fea- 
tures are typically used to measure the spectral types of 
M-type stars. This veiling, if not removed or if incom- 
pletely removed, will fill in the molecular absorption lines 
and lead to a systematic overestimation of the spectral 
types. 

By utilizing the [770 index] to derive T c g we are instead 
largely immune to these effects, and we are now able 
to correct the previous spectroscopically derived T c g via 
equation Equation [3j 

In conclusion, our technique to measure stellar T e g is 
able to derive a more accurate estimate of T e g than op- 
tical spectroscopy, being less biased by a possible veiling 
excess. 

3.4. Calibration of the [Ay index] and derivation of 
extinctions 



In section |3.1| we derived an empirical relation be- 
twee n th e [770 index] and the stellar T e g, and in sec- 
tion 3.2 we explored the effects of accretion luminosity 
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Fig. 11. — The color-color plot expressed in units of [770 index] (reddening i ndependent) and \Ay index] (parallel to the reddening 
direction). Left panel: sample of extinction-corrected sources; Ay values are from |Riddick et al.| ^2007^ , the JHK photometry of |Rqbberto| 
|et al.| | |2010[ ), and the optical analysis of Paper II. For the latter, accretors and candidate non accretors are distinguished as in Figiire|s] The 
solid line represents our best fit of the data, i.e. the empirical isochrone for Ay = 0. Right panel: same diagram, reporting all the sources 
in our catalog, with no correction for extinction. These are color coded according to their photometric errors. The entire population, with 
the exception of very few sparse sources, is located above the isochrone, at positive values of extinction. 




Fig. 12. — (/ — J) vs T c ff : the open dots are the measured colors, 
the filled gray dots are the dc-rcddcncd ones. Black dots are dered- 
dened colors for sources with know estimate of accretion excess, 
showing no significant ex cess. Green squares are sources whose 
T c ff is from Riddick ct al. (20071. The dashed line is the prediction 
from sy nthetic photometry on the BT-Settle atmosphere models of 
|Allard et al.||2010[ |; the solid line is the empirical fit to the data. 

on there derivations. Here, similarly, we find the relation 
between the [Ay index] , obtained in Equation [TJ and 
the actual extinction Ay of the ONC stars, in magni- 
tudes. In the color-color diagram, the extinction Ay of 
a star is the distance, along the reddening direction, be- 
tween the observed position and the isochrone, whereas 
the [Ay index] is the distance, also along the same di- 
rection and scaled by the appropriate A\/Ay, between 
the observed colors and the abscissa (770 — 753) = 0. 
Therefore the [Ay — index] and the true Ay differ only 
by a constant, which depends on the [770 index] (i.e., it 
is a function of T e g). 

To this end, we plot the extinction-corrected ONC 
stars wit h independent Ay es timates from Paper II 
and from Riddick et al. (2007) in the [Ay index] vs. 
[770 index] plane (Figure IT - left-hand panel). We also 



Robberto et al. (2010), from which we use stars with the 
lowest ^eff (low [770 index]) and with the smallest pho- 
tometric errors. For the co olest T s ff considered he re, the 
J vs (J—H) isochrone from |Robberto et al. ( 2010 ) is ver- 
tical, and its color (J — H ~ U.t>5) is independent of age, 
therefore Ay can be derived by simply dereddening the 
sources on the iso chrone. By not includi ng the i^-band 
photometry from Robberto et al. (2010), the infrared- 
derived extinctions should not be strongly affected by 
the possible presence of disk excess. The thick line in the 
figure is our fit to the data, i.e., the empirical isochrone 
corresponding to Ay=0. This also represents the offset 
to be subtracted from the measured [Ay index], as a 
function of the [770 index], to compute the true Ay. 

Among the points in Figure [IT] from our previous Pa- 
per II, we can distinguish accreting sources from non- 
accretors, just as in Figure|8] In Section |XT| we suggested 
that previous spectroscopically derived T c ff estimates are 
systematically overestimated for accreting members in 
the M-type range. This is now confirmed by Figure 
111! whereas the extinction-corrected non-accretors lie 
on a curve consistent with the lower edge of the non 
extinction-corrected population (right panel of the same 
Figure), for accretors we find Ay estimates systemat- 
ically overestimated. This is because if T e g is over- 
estimated, the intrinsic colors are underestimated, and 
therefore the color excess attributed to extinction is over- 
estimated leading to higher Ay. 

Figure |11| right-hand panel, shows again the empir- 
ically calibrated isochrone (i.e., for Ay = 0), together 
with all the stars in our catalog. Unlike in the left-hand 
panel, we have not subtracted the extinction here. Over- 
all the population lies above the line, i.e., at positive 
values of Ay. As would be expected if the isochrone is 
correct, the distribution of sources above the isochrone 
also appears uniform for different [770 index]. There is 
an overabundance of highly reddened stars at the high- 
est useful limit of the index ([770 index] = —0.1), i.e., 
at the highest limit of the T r range accessible with our 
medium-band photometry. This, as will be shown below, 



use the ISPI JHK photometry of the ONC presented in 
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4. THE H-R DIAGRAM 



Having derived two empirical relations, one to obtain 
the T e ff from the [770 index] for stars in the spectral 
subtype range M3-M8.5, the other to derive Ay once 
the [770 index] is known, we now derive the stellar pa- 
rameters for the sources in our catalog. 

Sources lying below the isochrone of Figure nib appear 
to have negative extinctions, which is unphysicaT; 93 stars 
showing an extinction in the range —2 mag< Ay <0 mag 
are therefore assigned Ay = 0, essentially shifting the 
measured colors to the closest physically meaningful so- 
lution. Another few stars with Ay < —2 mag were re- 
moved from the sample, being likely spurious detections 
or faint sources with very inaccurate photometry. The 
fraction of sources showing negative extinction (about 
5%) is much smaller than in previous works (e.g., Pa- 



per II or |Hillenbrand [1997 1 , validating the improved ac- 
curacy of our methods to derive stellar parameters. In 
total, we derived T e g and Ay for 1280 sources, whereas 
only 544 of them had a previously assigned T c g. Merg- 
ing the new T e g and Ay with those of our previous work, 
our sample counts 1807 sources with both T e g and Ay 
from which we can derive an updated H-R diagram for 
the ONC extending to well below the H-burning limit. 

4.1. Bolometric corrections 

In order to convert the de-reddened magnitudes into 
logL/L Q , knowledge of the bolometric corrections (BCs) 
as a function of T c g is required. Unfortunately, we can- 
not rely for such relation s on previous works (e.g., [Flower 
1996 Bessell et al.|1998| for several reasons: 1) the pho- 
tometric bands used here, including /—band, are non- 
standard; 2) optical BCs usually do not extend in the 
BD temperature range; and 3) if the optical stellar col- 
ors are age dependent, as suggested in Paper II, the BCs 
may be also. Thus, BCs valid for main-sequence dwarfs 
might not be adequate for a young populations such as 
the ONC. Instead, we use synthetic photometry to derive 
our BCs, according to the formula: 



BC = 2.5 • log 



Jfx(T)-S x d\ / J f x (Q) • S x d\ 
Jfx(T)-d\ I Jfx(&)-dX 



, (4) 



where f\(T) are synthetic spectra, and Sx the filter 
throughput. 

As we have shown above (see also the discussion in 
Paper II), current atmosphere models do not correctly 
predict the optical broad band colors for cool stars. This 
clearly affects the computation of BCs, in particular BCs 
in /—band, BC\. On the other hand, the bolometric flux 
for a given synthetic spectrum will be trivially correct at 
a given T e g since, by definition, L oc R 2 T A . The radii 
R, which could also be inaccurate from the evolutionary 
models, do not affect the BCs since they introduce an 
identical proportionality factor in both the numerator 
and denominator of equation [4j 

From the comparison of our data with the synthetic 
isochrones (see Figure [6]) we cannot ascertain the cor- 
rectness of the synthetic /-band fluxes. It is true that 
the isochrones are unable to reproduce the data in the 
color-color diagram, but the shortcomings may be princi- 
pally in our medium bands and not in the /-band. At the 



same time, the BT-Settle models ( Allard et al. 2010 ) have 
been validated in the NIR for dwarf stars, such that the 
predicted J and K fluxes match the observations down to 
BD masses. Thus, if the predicted /—band magnitudes 
are also correct, we would expect the predicted (/ — J) 
color to be also consistent with the observations across 
all T eff . 

We perform this test using the J-b and photometry of 
the ONC from Robberto et al. (2010), toge ther with our 
/—band magnitudes and (as in Section 3.4) the Ay esti- 
mates from previous works. Figure fl2] shows the dered- 
dened (/ — J) color as a function of T c g for our ONC 
sample, in comparison with the same quantity computed 
usi ng the BT-Settle atmospheres, assuming a Baraffe et 
al. (19981 1 Myr isochrone to constrain the surface grav 



lty as a function of T e g . As in Figure [T2j we highlight 
sources showing no or little accretion excess, except at 
the coolest T e g (T e g < 3000 K), for which we do not have 
estimates of accretion from Paper II. For T e g > 3200 K, 
the isochrone fits the data well. However, the match is 
worse at lower T e g, the isochrone predicting lower (I— J) 
than shown by our data. We checked that this conclusion 
is not strongly dependent on log g. 

Trusting the predicted J— band fluxes, as explained 
above, we conclude that the BT-Settle atmosphere mod- 
els tend to systematically overestimate the /—band 
fluxes, and this effect is larger for cooler T e g. In Fig- 
ure [l2] we also show the empirical fit to the data. The 
difference between the two lines represents the empirical 
correction needed to calibrate the predicted /-band syn- 
thetic photometry using the BT-Settle models. Whereas 
this is small for T c g > 3000 K, in the BD T c g range the 
correction increases up to 1 magnitude. Thus, we correct 
the computed BCs, derived from Equation |4j applying 
this offset as a function of T e g. 

Finally, we derive the bolometric luminosities for all 
our sources from: 

log£/£ =0A-(M I<e -I+A I + DM + BCi(T eS )) (5) 

where Mr )0 = 4.03 is the absolute magnitude, in WFI 
/ band, of the Sun (Paper II) and DM=8.085 is the dis- 
tan ce mo dulus, adopting d = 414 ± 7 pc from Mentenet 



al. A2007I 



le resulting H-R diagram is shown in Figure |13| 
where we show separately the stars whose parameters 
have been derived from the new data presented in this 
work, and those, generally with T c g > 3500 K, whose 
T e ff and L are from Paper II. The apparent discontinu- 
ity in the number of sources around logT c ff sa 3.54 is 
mainly due to the much higher completeness in the num- 
ber of low-mass stars with T c g estimates provided by the 
present work (see Section 4.3). 



4.2. Background contamination 

Despite the large column density of the Orion 
molecular cloud, which produces an extinction up to 
Ay = 100 mag (iBergin et al. 1996 Scandariato et al. 



2011 1, contamination from background stellar popula- 



tions starts to become non negligible at faint magnitudes 
in t he ONC. In fact, as suggested by previous works. 



e.g.lRobberto et al |20l6) [Hillenbrand fc Carpenter|2000 
Muench et al.|2002[) , the traction ot background contain- 
er sus membi 



mants versus 



)ers increases towards the substellar 
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Fig. 13. — The new H-R diagram of the ONC. Green dots are the new stars whose stellar parameters have been derived in this work, 
blue dots are members from Paper II, mostly having a large enough T B g to not be derivable using the [770 index], or saturated in our new 
WFI photometry. 
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Fig. 14. — Stellar density map in the H-R diagram. The dashed 
line marks our adopted separation between bona fide members and 
contaminants. 

mass range in the ONC. This is both due to the increase 
in the number of faint background sources seen through 
the Orion molecular cloud, and the decrease of actual 
members of the young cluster in the BD mass range, 
because of the declining initial mass function (IMF) at 
substellar masses (see below). 

The contamination from background sources is evident 
as a large number (>200 stars) of faint sources in our H- 
R diagram (Figure 13 1, which arc well separated from the 
cluster sequence. Given the relatively low T e g for these 
objects (~ 3500 K< T cff <~ 2900 K), they could be 
either cool dwarf stars or red giants, the latter more eas- 
ily detectable at large distances and through the Orion 
Molecular Cloud. The relative number of these sources, 
relative to young cluster members, becomes dominant at 
our detection limit. 

A small number of these relatively less-luminous 
sources might be actual ONC members, whose flux 
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Fig. 15. — Spatial distribution of the ONC members placed in 
the H-R diagram (black dots). The red circles mark the position 
of the faint sources identified as background sources. 



may be masked (e.g., by edge-on circumstellar disks), 
and therefore are visible only in scattered light. This 
phenomenon has been o bserved in several young star- 
forming regions (s ee e.g., Guarcello et al.|[2010 Kraus & 



Hillenbrand 2009 



and references therein); however the 
expected fraction of such objects is very smal l, including 
no more than a few percent of the population (|De Marchi 



et al. 2011). On the other hand, the faint sources we 
detect below the young sequence are about 16% of the 
total number of stars with T c g < 3200 K. Therefore, we 
can assume that most of these sources are actually con- 
taminants. Moreover, removing a small fraction of faint 
members erroneously considered as background contami- 
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Fi g. 16. — The H-R dia gram of the ONC members, after the removal of contamination from background stars. Evolutionary models are 
from [Baraffe et al.| l |1998l l (upper panel) and |D'Antona fc Mazzit clli (1998 ) (lower panel). 



nants will not affect significantly the IMF derived in this 
work, given the small expected number of these objects 
compared to the entire studied population. 

We proceed to attempt to remove the background con- 
tamination from the H-R diagram. To this purpose, we 
define a cut in the H-R diagram such that stars satisfying 
the following criterion are considered as non-members: 



tion provided by the Orion Nebula is hig her on a central 



no rth-south strip, centered o n the ONC ( Scandariato et 
[al] |2011| |Bergin et ah] |1996[ ) . Therefore the density of 



logL < -1.93 

logL < -5.607 

+65.848 -(logT cff 
-253.02- (log T cff - 



logT cff > 3.54 



- 3.4) 

3.4) 2 logT t 



off 



< 3.54 



This is shown in Figure [14] All sources located below 
this line arc regarded as non -members and excluded from 
further analysis. Figure 15 presents the spatial distribu- 
tion of contaminants in comparison with that of the can- 
didate members. Whereas the latter are centrally con- 
centrated, with an evident north-south elongation, con- 
taminants appear uniformly distributed. Also, we find a 
slight overabundance of contaminants on the east side of 
our field of view. The reason for this is that the extinc- 



Background objects is expected to be higher on the two 
sides of this band. Here, the nebular emission of the 
OMC is significantly brighter on the west side of our 
FOV, and this decreases the detection limit in this area. 
Therefore the east edge of our FOV is expected to show 
the most prominent density of background contaminants. 

In Figure [16] we present once again the H-R diagram 
of the ONC, now including only our candidate members, 
together with PMS evolutionary models. In particular 
we use is ochrones and tracks from eith er Baraffc ~et al.l 
( 1998 ) or ID'Antona fcM azzitcllil ( |1998[ ) These models 
unlike e.g., the (5TSs = eT 5 aT hf ' ' ! n 11 - 1 Cl 



W\) OTThe |Palla fc Stahler 
( |1999| , also extend into the substellar mass range, and 
therefore are more suitable for our sample, which reaches 
masses as low as 0.02 M Q (20Mj up ). 

4.3. Completeness in the H-R Diagram 



In Section 2.3 we have derived the photometric com- 



pleteness of our survey. Starting from this result, we now 
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Fig. 17. — Completeness function in the H-R diagram, computed accounting for photometric detection and differential reddening as 
described in the text. The left panel shows the result as computed with our monte-carlo simulation, valid for the spectral types obtained 
from the medium-band photometry analyzed in this work, for logT c ff < 3.54 (T ff < 3500 K). The right panel shows the correction of this 
completeness after accounting for the additional stars placed in the HRD whose stellar parameters are taken from Paper II. 



evaluate the completeness function in the H-R diagram. 
To do this, we transform the H-R diagram (i.e., T ff and 
logL) into our observational space (i.e., the 753, 770, and 
/-band magnitudes) , and assess the extent to which red- 
dening causes sources of various intrinsic T e g and log L 
to be missed in our survey. 

We set a uniform, dense grid of logT e ff and logi bins 
in the H-R diagram, limited to the region logT G ff < 3.54. 
For each bin, we apply in reverse the relations derived 
in Section |3.1| and |3.4| to derive from T e g and log L the 
[770 index], the [Ay index], and (using the BC) the 
unreddened /-band magnitude, Iq. Then, by inverting 
equation [2] and [T] we also derive the magnitudes 771753 
and 777,770. Thus we have associated to each pair of H-R 
diagram parameters (logTeff, log L) the triplet of intrinsic 
magnitudes (777753, 777770, /o) f° r Ay = 0. 

Next we apply reddening to these magnitudes using the 
empirical reddening distribution of the ONC stellar pop- 
ulation. As in Paper II, we limit ourselves to the most lu- 
minous half of the ONC population, as these stars sample 
the full range of Ay. With a Monte Carlo approach, we 
draw 1000 Ay values from the reddening distribution and 
apply them to the intrinsic magnitudes (777753,777770 ,/o) 
associated to each point of the H-R diagram. For every 
i-th (1 < i < 1000) triplet of reddened magnitudes, we 
compute the photometric completeness (7753(1), Crro(i), 
Ci(i). The product of these three, C(i), averaged over 
the 1000 values of Ay , provides the "differential redden- 
ing normalized" completeness for that particular bin in 
the H-R diagram. 

The result is shown in Figure 17 1. We find that the 



completeness is mainly dependent on logi, while its de- 
pendence on T b q for a given luminosity is small and ap- 
pears only at low temperatures (logT^ < 3.45). We 
note also that the completeness is rather high (~ 50%) 
at the low-mass end of our catalog, and thus even with- 
out correcting for the completeness, the observed number 
of ONC members in the substellar mass range is correct 
to within a factor of ~2. 

Our H-R diagram also includes sources whose stel- 
lar parameters are taken from Paper II instead of our 
new medium-band photometry. For all members with 
logTcff' > 3.54, we simply adopt the completeness func- 
tion in the HRD derived in that paper. In addition, there 



is a group of stars with logT c ff < 3.54 but which are not 
present in our new medium-band photometry and thus 
whose stellar parameters are from Paper II. The presence 
of these sources increases the overall completeness in this 
region of the HRD relative to that computed as above. 
To account for this, we divide the HRD into a grid with 
spacing 0.01 dex in logT e ff and 0.1 dex in logL, and we 
count in each bin the relative number of sources from 
Paper II versus members classified from our new data. 
This ratio, smoothed in logT e ff and logi, is added to 
the completeness shown in Figure [PT^ . 

The final result is shown in Figure |17b . The evident 
discontinuity in completeness at logT c ff = 3.54 is due to 
the fact that for larger T c ff the stellar parameters are ob- 
tained from optical spectroscopy on an incomplete frac- 
tion of sources. From the two-dimensional completeness 
shown in Figure |17b , we are able to assign a complete- 
ness correction to each our sources, allowing us finally to 
derive a completeness corrected initial mass function for 
the ONC. 

5. THE INITIAL MASS FUNCTION 

To derive the initial mass function (IMF) of the ONC 
from our observations, we use pre-main-sequence evolu- 
tionary models to convert our T e g and log L from the 
H-R diag ram into masses and ages. W e do this using 



al 
mo 



both the D'Antona & Mazzitelli (1998) and Baraffc ct 



1998 i models (Figure 



els, for a large fract 



re |i"6| 
ion - (a 



In the case ot Baratte 
about 25%) of sources we 



cannot assign masses and ages, as these stars are located 
above the 1 Myr isochrone, the minimum age computed 
for this family of models. This not only decreases our 
stellar sample, but also biases our findings, in partic- 
ular the mass distribution. To overcome this selection 
effect, we consider two cases for the mass estimates from 
Baraffe models: a) we reject all source above the 1 Myr 
isochrone; b) we include them by assigning a mass based 
on the Teff-mass relation of the 1 Myr isochrone. Since 
for very-low mass stars the PMS evolutionary tracks are 
nearly vertical in the H-R diagram, this approximation 
is fairly good. In Table [4] we present the derived stel- 
lar parameters for the ONC sources, using both sets of 
models. 

From these masses, we derive the mass distribution 
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Fig. 18. — The measured IMF for the ONC, fitted with a log-normal distr ibution (left panel) o r a t wo-phase power law (right pan el). The 
top and bottom panels represent the mass distributions obtained assuming |Baraffe et al.| | |1998| and |D'Antona fc Mazzitellij l |1998| l models 
respectively, whereas the center panels show the result assuming Baraffe models and including also stars located above tne 1 Myr, whose 
mass has been extrapolate d from the ir T e fj (see text ). The shaded areas enclose the 90% confidence interv al for each fit. T he red curves 
represent the IMF of |Chabrier] | |2003| (left panel) and |Kroupa| j2001| (right panel); the red dashed line is the[Chabricr (2003|l system IMF. 



TABLE 3 
IMF 





lognormal 
logrfi c cr(logm) 


2-phase powerlaw 
Ti T 2 logm c 


BCAH98 
BCAH98 all stars 
DM98 


-0.45 ±0.02 0.44 ±0.05 
-0.53 ±0.02 0.39 ±0.03 
-0.56 ±0.02 0.38 ±0.01 


-1.12 ±0.90 0.60 ±0.33 -0.53 ± 0.26 
-3.05 ±0.53 0.39 ±0.10 -1.00 ±0.08 
-2.41 ±0.25 1.30 ±0.09 -0.77 ± 0.02 



Note. — The power law exponents Y follow the standard for which the Salpeter slope is 
r = 1.35 



£ (log m) by binning the ONC members in equally-spaced 
mass bins. For each mass bin, we account for its exact 
completeness by adding the inverse of the completeness 
of each source. We associate an uncertainty distribution 
to each measured value of £(logm) equal to a Possion 
distribution of mean fi — Ni, where Ni is the number of 
sources in the i— th bin, scaled by a factor equal to the 
overall completeness correction for that bin. We stress 
that, strictly speaking, our mass function is actually a 
"system" mass function rather than a proper "initial" 
mass function, in the sense that we do not account for 
unresolved binaries or multiple systems. This, however, 
does not influence significantly our results, since the bi- 



binaries) in Orion is small (< 15%, Padgett et al. 


1997 


Petr et al. 


1998 


Reipurth et al. 


2007 j, and about half 



in our observations. 

It it well established (e.g., Bastian et al.|2010 ) that the 
IMF generally follows a power-law in the intermediatc- 
and high-mass range (M > M ), whereas for low-mass 
stars and BDs — which is the region of the mass spectrum 
most relevant for our study — this function can be approx- 
imate d with a shallo wer power law £(logm) oc m _ ^ r+1 - ) 
(e.g.. ' 



£(logm) 



Kroupa||2001|) or with a log-nor mal distribution 

log rn 



-(log? 



72<7 a 



(e. 



Chabricr 



We use both forms to fit our measured IMP in the 



2003). 
TJNTJ. 



We use a M onte Carlo simulation following Da Rio 
et al. (2009a) to account for the uncertainties in the 



measured star counts, as follows. For every mass bin, 
we consider the error bars with their statistical distri- 
bution, and generate a large number [n — 10, 000) of 
points drawn from the error distribution. Then an un- 
weighted fit is run on all these (n times the number of 
bins) points. The best-fit parameters are isolated us- 
ing a Levenberg-Marquard minimization algorithm. The 
uncertainty associated to the parameters has been com- 
puted with a sampling technique as follows: for every one 
of the m bins, we randomly consider only one of the n 
values previously simulated, i.e. a random sample from 
the distribution describing the error bar of the bin, and 
we fit the IMF function on these m data points, deriving 
the best-fit parameters. By iterating this selection and 
fit process 1000 times, we derive 1000 sets of parameters. 
The standard deviation of each parameter for the 1000 
tests is the uncertainty in the estimate of the parameter 
itself. 

The resulting fitted functions for both sets of evolution- 
ary models are shown in Figure [18] In Table [3] the fitting 
parameters are reported: the characteristic mass and the 
width for the log-normal fit, and the two power law slopes 
as well as the break point for the 2-phase power law. We 
find that the two families of evolutionary models lead 
to significant differences in the derived IMF. Whereas 
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the Baraffe tracks produce a smooth distribution, which 
appears well fitted by a log-normal distribution charac- 
terized by a continuou s change in the IMF slope, the 
|D'Antoiia &: Mazzitelli models lead to a more peaked dis- 
tribution, with an evident maximum at M ~ 0.15 M . 
The reason for this difference can be attributed to the 
different shapes of the evolutionary tracks in the very 
low mass range (0.3 M > M > 0.1 M ), where the 



models of |D'Antona fc Mazzitelli predict a larger area 
of the HRD covered by adjacent track s. The shape of 
the IMF derived from the IBaraffe et al.l models does not 
depend significantly on which method we adopt for in- 
clusion of very young stars lying above 1 Myr isochrone 
(see above). 



In Figure 18 we also compare our measured mass dis- 
tribution with standard refere nce IMFs . In p articular 
we consider the disk IMF from Chabrier ( |2003[ ) which is 
expressed as a log- normal function, with a characteristic 
mass of either 0.079 M Q or 0.22 M respectively for the 
"single objects" IMF or the "system" IMF, the latter not 
corrected for stel lar mult i plicity . Also, we consider the 
universal IMF of Kroupa (2001 1, which is a multi-phase 
power law with shallower slopes in the VLMS and BD 
regime. Depending on the assumed evolutionary mod- 
els and the type of function fitted, our IMF peaks in 
the range 0.1 M — 0.3 M , the higher values being 
larger than that most often reported, even in the case of 
unresolved multiplicity. This finding, as well as the high 
dependence on the evolutionary model family used to de- 
rive stellar masses, is consistent with what we have found 
in our previous work (Paper I) limited to the stellar mass 
range. 

A striking feature of our IMF is the steep decline in 
the BD mass range (see Figure 18); here our mass distri- 
bution is much lower than the standard IMFs, implying 
a significant deficiency of substellar objects in the ONC. 
We stress that this finding is not biased by detection in- 
completeness, given that the derived mass distribution 
has been corrected for completeness star by star, and 
the completeness function we have derived in Section 
4.3 remains fairly high (~ 50%) even in the BD mass 
range. Moreover, even hypothesizing an overestimated 
completeness for the smallest masses in our sample, the 
relative scarcity of substellar members is evident up to 
about 0.1 M ; there simply are not many very cool, very 
low luminosity sources observed. Such a steep decline of 
the mass function in the BD mass range differs also from 
IMF determinations in other young regions. In fact the 
majority of young clusters have been reported to show 
substellar mass function compatible with the Kroupa or 
the Chabrier mass dis t ribut ions. Some examples are 
Chamaeleon I |Luhman| ( |2007[), Uppe r Sco (Lodieu et al 
20071), the a Ori region lUaballero||2009|), A ( Jri ( |Bayo 
eTlilpmj ), NGC 6611 floiiveira et al.||2009|). SmufaF 



findi ngs have also been r eported tor oider open clusters 
(e.g ., iBouvier et al.||2005 ); in fact for the Praesepe clus- 
ter, | Wang et al. 1 201 1 ) measure even an increasing IMF 
down to ~ 70 JVl j . 

Previous studies of the IMF in the ONC also found a 
decrease in the relative number of members in the sub- 
stellar mass range, both in absolute number and rela- 
tive to the s tandard IMFs. Studies based on photometry 



alo ne (e .g., Hillenbrand & Carpenter 2000 Muench et 



aT]|2002 ) using the conversion of near infrared LFs, find 
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a measured IMF slope below the H-burning limit not 
steeper than T ~ —1, which is the shallowest value we 



derive assuming Baraffe et al. ( 1998 1 models, and up to 2 
units natter than what we derive with D'Anto na &: Mazz^l 
itellil A 19981) tracks. 



Also, 



Slesnick et aL] ( |2004[ ) finds an 



1Mb for the ONC rapidly decreasing below the H-burning ONC of about l(r M Q . 
li mit, but then flatteni ng in the BD mass range. 



works, which we believe is much improved through the 
methodology that we have developed here. 

Based on our measured IMF, complemented in the 
intermediate- and high mass range using the results from 
Hillenbrand (19971, we assess a total stellar mass for the 



Muench et al. ( 2002 ) reports a secondary peak in the 
mass distributi on at about the deute rium burning limit 
(~ 0.17 M ); |Slesnick et"al] p004| tentatively find a 



similar feature, although with less significance and at a 
slightly higher mass. This mass roughly corresponds to 
our detection limit. However, from our results, we do 
not find any evidence of either a flattening of the IMF 
or an increasing number of sources for the lowest mass 
bins. We interpret this inconsistency as an inaccurate 
estimate of the background contamination in previous 
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